This example is based on the simple model of two intersecting faults used by Mark in some of the presentations. Idea: we know that two faults exists, but not their age relationship. Offset is expected/ potentially estimated from stratigraphic relationships across fault.
In [17]:
import sys, os
import subprocess
noddyprogram = os.path.realpath(r'../../noddy')
example_directory = os.path.realpath(r'../../docs/examples')
# set Pythonpath to pynoddy
sys.path.append(os.path.realpath(r"../../python"))
import pynoddy
In [9]:
# Compute noddy model for history file
history_file = 'fault_model_1.his'
history = os.path.join(example_directory, history_file)
output_name = 'noddy_out'
# call Noddy
print subprocess.Popen([noddyprogram, history, output_name],
shell=False, stderr=subprocess.PIPE,
stdout=subprocess.PIPE).stdout.read()
#
In [15]:
In [16]:
In [ ]: